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We study the interplay between graviton reggeization and double- 
logarithmic in energy contributions to four-graviton scattering in theories with 
and without supersymmetry. Predictions to all orders in the gravitational cou- 
pling are given for these double-logarithms. As the number of supersymmetries 
grows these terms generate a convergent behaviour for the amplitudes at very 
high energies. At two-loop level, we find agreement with previous exact results 
for Af = 8 supergravity and with those of Boucher- Veronneau and Dixon, who 
studied the Af = 4, 5, 6 supergravities using a conjectured double-copy struc- 
ture of gravity. 

1 Introduction 

Collider phenomenology and the anti de Sitter / conformal field theory 
(AdS/CFT) correspondence [l}j3] have motivated numerous investigations of 
scattering amplitudes in gauge and gravitational theories in recent years. In 
this context, the study of the Af = 4 supersymmetric Yang-Mills theory has 
been particularly successful since it is possible to investigate loop corrections by 
only evaluating a reduced set of master topologies |4j. These results for Af = 4 
SUSY can then be used to also obtain amplitudes in Af = 8 supergravity, 



1 



where the set of minimal topologies is still valid, and to address the question 



. This procedure is 
which has recently 
In Einstein-Hilbert 



of the renormalizability of the theory at higher orders [5j 6 
based on a conjectured double-copy structure of gravity 
been applied to J\f = 4, 5, 6 supergravities at two-loops 
gravity progress is slower (9-11 since supersymmetry or string theory based 
techniques 12 14 cannot help. 

It is also possible to get important information of the all-orders structure of 
scattering amplitudes when they are considered in certain kinematical regions. 
An interesting example is the study of graviton scattering in multi-Regge kine- 
matics (MRK). In this case the amplitudes present a factorized form which 
can be interpreted in terms of the exchange of reggeized gravitons 15 16 
together with eikonal and double-logarithmic terms 17 19 . These contribu- 



tions, together with new interaction vertices, can be described by means of a 
high energy effective action p20j. It is noteworthy that the graviton emission 



vertex in MRK can be written as a double copy of the corresponding 21 23 



QCD gluon emission vertex 17 19,24 



In the present work double-logarithmic in energy contributions to four- 
graviton scattering to all orders in the gravitational coupling will be evaluated. 
This will be done for arbitrary supergravities as well as for Einstein-Hilbert 
gravity. We will improve previous results based on the resummation of ladder- 
like diagrams by considering the full set of contributing topologies. The trun- 
cation of our results to two loops is in agreement with recent calculations in 
the literature for TV" = 4, 5, 6, 8 supergravities. The all-orders resummation of 
these contributions generates amplitudes which grow with energy when < 4 
and asymptotically tend to zero when M > 4. M = 4 supergravity corresponds 
to a critical theory where the leading double-logarithmic contributions cancel. 



2 Double- logarithmic approximation 

For our analysis it is convenient to define the following normalization for the 
four-point amplitudes 

Ai,(N) = A^t^M^N), (1) 



A*™ = k 2 ^, (2) 
tu 

oo 

M m) = 1 + J2 M W)> ( 3 ) 

L=l 

where L corresponds to the loop order and N labels the number of gravitinos 
in the theory. The Mandelstam invariants are s — {jp\ + P2) 2 , t = (p\ — ps) 2 
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and u — (pi — P4) 2 - k 2 = SttG, with G being the Newton's constant. 



2.1 One- loop amplitudes 

When the one loop amplitude -M^ 1 ^ is calculated in the Regge limit (with 
s ^> — t = \q\ 2 ) the graviton Regge trajectory 17-19 



«M 2 f M I A q- kf + (k q- kf _ |g|2 + N ft _,_ k) 



vr J \k\ 2 \q- k\ 2 \ \k\ 2 \q-k\ 2 ™ 2 

appears multiplied by ln(s/|g| 2 ). We have used the notation a = k 2 /(8ti 2 ). 
This expression contains both infrared and ultraviolet divergencies which can 
be regulated by, respectively, the cut-offs A and A, to obtain 



VI 



2 N - 4 A 2 



= -a\q\ 2 ( v ln ^ + 11 y^ ln Jqf) ' (5) 

The ultraviolet divergence at A — > oo is not fundamental because gravity is 
renormalizable at one loop. It has a kinematical origin which means that the 
parameter A should be substituted by y/s, leading to the appearance of the 
double-logarithmic term ~ a In 2 s in the elastic scattering amplitude. 

To understand this point in more detail, let us recall that the one loop 
contribution at high energies can be obtained making use of the Sudakov 
parametrization for the virtual particle momentum, 

I ^ I 

k = (3 pi + a p 2 + k± , d 4 k = — dad(3 d 2 k± , (6) 

where p\ and p 2 are the momenta of the colliding particles. Calculating the 
Feynman integral over a by residues we obtain with leading logarithmic accu- 
racy the following expression 

(1) If d 2 k ± f s d(J3s) 



^ M)~- / TTrf 2 / o ,, V (7) 

J \k ± \ 2 J\ q \2 (3s - IfcJ 2 

It is then clear that the ultraviolet divergence at \k±\ — > oo is absent. In 
the infrared region of integration \kj_\ \q\ 2 the above expression factorizes 
in the form ^l^ 1 '' ~ ln(|g| 2 /A 2 ) ln(s/|g| 2 ). This Regge factor containing the 
infrared divergence can be extracted to all orders in perturbation theory and 
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the amplitude with double-logarithmic accuracy can be conveniently presented 
using the following Mellin transform: 



where the tree level amplitude is A® orn ~ n 2 s 2 /\q\ 2 and the integral over u 
contains the contribution of the virtual gluons and gluinos only with \kj_\ 2 > 
\q\ 2 . The t-channel partial wave in the double- logarithmic approximation 
can be expanded order by order in perturbation theory, 

A m = '£c^ (A) , ( 9) 

n=0 ^ ' 

where b is the dimensionless parameter 

b = a\q\ 2 . (10) 



In Ref. 17 one of the authors of this work (L.N.L.) made the assump- 
tion that double-logarithmic contributions in gravity appear only from ladder 
diagrams, which allowed him to obtain a closed expression for the scattering 
amplitude in terms of a Bessel function. However, even in simpler field theories 
like QED and QCD there is another source of double-logarithmic terms. For 
example, in e + e~ backward scattering the diagrams with virtual soft photons 
emitted and absorbed by the external fermions are essential. These contri- 
butions contain, apart from the integration region \k±\ 2 > \q\ 2 , the universal 
infrared divergencies from the region \kj_\ 2 <C \q\ 2 in the form of Eq.Q. It is 
then natural to realize that similar non-ladder contributions also play a role in 
gravity. At one loop they would lead to the following correction to the elastic 
amplitude 

41ft , f s d\k ± \ 2 f s d{(5s) 



- b f 


d\k ± \ 2 


J X 2 


\k ± \ 2 


— b In - 







^ orn \k±\ 2 J W (3s-\k^ 

which, together with the corresponding ladder correction (including gravitino 
contributions) , 

"gladder i. f ^ ~ ^ \ , 2 S 



A^ mn ~ I 4 / \q\ 2 
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generate the complete one loop correction (cf. Eq. (|5 

^-Hn^^^ln^,. (13) 



Al om \q\ 2 V A 2 V 4 y |g| 2 

To better understand this result let us now compare it to the exact one- 
loop amplitudes available in the literature. In Ref. [8] Boucher- Veronneau and 
Dixon used a conjectured double-copy structure of gravity to evaluate four- 
point scattering amplitudes at two loops in Af = 4, 5, 6 supergravities. They 
also calculated again the well-known Af = 8 case. 

For Af = 8 it is possible to write the exact one-loop amplitude (using 
s + t + u = 0) in the form 



M { }}„-* = a tin I — I In 



4,(iv=8) - V — t / V —t 



—s \ . / —u 



Double Logs 



Trajectory 



Eikonal 

We have indicated the terms which will generate, in the Regge limit (s ~ —u), 
the double-logarithms, the graviton trajectory and the eikonal pieces. 

Following Ref. [8] it is possible to relate the previous expression to the 
Af = 4, 5, 6 supergravity amplitudes, i.e. 



S / V S / V S J V s 



M$ N=5) = M$ N=8)+ at^(y(^y^), (16) 

m$ n=6) = M:u 8) +«4K ln2 (^) +7r2) - (17) 

From the work of Dunbar and Norridge in Ref. [12] we also know the one-loop 
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amplitude in Einstein-Hilbert gravity (N = 0): 



with 

Q{x) 



+ 



+ 



M 



(i) 

4,(iV=0) 



ili 1 u _ 



(4 - lOx + 2x 2 + 15x 3 - 5x 4 - 3x 5 + x 6 ) 



n-1 



x i ln z fx) + 2a;ln(a;) + 7r 2 + ^a; n [ ^ In (g) + ^ ^ ^ 



n=2 



341 437 47 2 37 3 r 4 n 5 x # , . , ^ 



30 30 X 2 X 



n=l 



X 

n 



961 97 85 2 o 4 

1 x x - 2x + x 

90 12 12 



(19) 



Using these expressions we can see that the double-logarithmic contributions 
to these amplitudes can be written in the form 



M 



(1),DL 
4,(JV) 



N-A\ fat 
2~ 



In 2 



s 

~t 



(20) 



in agreement with our result in Eq. (13). 



2.2 All-loop amplitudes 

Since the sources for the double-logarithmic contributions in gravity are the 



same as in gauge theories (see Ref. 25 27 for related calculations in QED 



and QCD), it is possible to write a similar infrared evolution equation for the 
partial wave introduced in Eq. pi), namely, 



(TV) 



a fW /N-6\ f (Ny 



which, in pictorial terms, stems from the equation 

X>XX 



(21) 




.(22) 



At the right hand side of Eq. (21) the first term proportional to b describes 
the contribution of the virtual graviton with the smallest value of k±. The 
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second (ladder) term describes the contribution from the pair with the two 
softest gravitons or gravitinos exchanged in the t-channel. Let us indicate that, 
generally, the emission of the virtual soft graviton with transverse momentum 
k± changes the momentum transfer (q — > q—k) for the basic scattering process, 
and hence modifies the power factors \q\ 2n — > \q — k\ 2n of the corresponding 
amplitude in each order of perturbation theory. Nevertheless, it is correct 
to neglect the corrections ~ kj_ because these terms cancel the logarithmic 
contribution appearing from the integration over kj_. 



It is important to indicate that the coefficients of Eq. (21) in front of the 



two terms proportional to b are chosen in such a way as to reproduce the 



one loop contribution calculated in Eq. (13). The perturbative solution of the 



infrared evolution equation (21) has the form 



{N) _ b(N - 4) b\N-4)(N-3) 

b 3 (N -A) (5iV 2 - 2QN + 36) 



+ 



8w° 

b 4 (N - 4) (7iV 3 - 47iV 2 + 118iV - 108) 
8w* 

b 5 (N - 4) (21iV 4 - 16(W 3 + 556iV 2 - 960iV + 648) 



I6w 
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+ ... (23) 



which leads to the following double-logarithmic asymptotics of the elastic am- 
plitude: 

A UN) (s,t) = A Bom ^J ® (N) (0, (24) 

where £ = a \t\ In 2 -Aj and 
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<«> K ) = x _(iV-4) £ (iV-4) e 
KJ 2 2 2 v ; 4! 

(A ~ 4 )- (5^-26^ + 36)^ 

6! 

+ {iy ~ V (7N 3 - 47iV 2 + 118iV - 108) 



(N 


S 

-4) 


(N 


S 

-4) 



e 

8! 

t (21iV 4 - 160iV 3 + 556iV 2 - 960iV + 648) + . .(25) 

Before investigating in detail how the amplitudes behave at the light of this 
all-orders result, let us compare its two-loop truncation with other calculations 
in the literature. 
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2.2.1 Two-loop truncation and comparison with iV = 4, 5, 6, 8 results 

In the notation of Eq. ([3]), it is well-known that the complete two- loop ampli- 
tudes can be written as the sum of two pieces: 



m 



(2),N 



+ 



(26) 



where m^' N is infrared finite and has been calculated by Boucher- Veronneau 
and Dixon in Ref. Isl for N = 4,5,6. The N = 8 case was calculated earlier 



in Ref. 12 . We will now present these infrared finite remainders in the Regge 



limit and with leading double- logarithmic accuracy. First, it is convenient to 
write the exact amplitudes in the form 



A 



4, (AT) 



A 



Born 



X 



1 + 



-t 
A* 

N ■ 



a t 

T 

2 



iV-4 



In 2 

a t 

T 



In 4 



s 

~t 



m 



(2),N 
4,DL 



+ 



(27) 



where we have exponentiated the infrared divergent terms and singled out 
the non-exponentiating double-logarithmic contributions. The latter contain 

.ast 
il: 



those pieces related to the one-loop result and its square, following the 

J 4,DL 
2 



term of Eq. (|26|), and the mf]^ contribution, which can be read off Ref. 



m 



(2),AT=4 
4,DL 







m 



(2),AT=5 
4,DL 



m 



[2),N=6 
4,DL 



m 



(2),N=8 
4,DL 



at 

1 fat 
24 \ Y 

at 

T 

1 fat 

3 It 



In 4 



In 4 



-t 

s 

~t 







In 4 



In 4 



-t 



0. 

6 
1 

2 
5 
3 



a t 
T 

at 

T 

at 

T 



In 4 



In 4 



In 4 



s 

~t 



-t 



(28) 
(29) 
(30) 
(31) 



At the right hand side of these expressions we have written the final double- 
logarithmic contribution to the amplitude. It is important to note that these 
results are in complete agreement with the two-loop truncation of our predic- 



tion for any iV in Eq. (25). Higher order terms for different supergravities or 
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Einstein-Hilbert gravity can be obtained from it. This should serve as a useful 
test of multi-loop calculations of four-graviton amplitudes. As an example, in 
Af = 8 SUGRA we obtain 



°*0»te)+**(f)) 

,iBorn I " \ 

4,(JV=8) 



A.AJM = S\ — ^4?° m 



2 

at \ , 9 / s \ 5 / a £ \ , A / 5 



37 / a 6 f s \ 353 / a A 4 j 8 / s 

Now we turn to study the high energy asymptotic behaviour of the resummed 
amplitudes. 

2.2.2 Resummed supergravity amplitudes at high energies 

In terms of double-logarithmic contributions the simplest amplitude is that of 
Af = 4 SUGRA since their contribution adds to zero and we have the pure 
Regge asymptotic behaviour 



Aux-u(--t)= Af m (^) • (33) 



-a\q\ 2 lnlf£ 



In the case of M = 6 SUGRA Eq. (21 ) can be solved in the form 



OO 9 
z 2 b 



fi N=6) = dze-*e-i^, (34) 



obtaining the following result for the amplitude 



I 1^ 

A,, (N=6) (s,t)=AT in ^ ^exp^-^ln 2 ^. (35) 

In the general case with arbitrary iV it is useful to introduce the new function 
y(x) and the new variable x according to the definitions 

/«">— = (36) 
9 



to reduce our Eq. (21) to the Riccati equation 



yW\ X ) + yW\ X ) - X yW + = . (37) 



By introducing the new function tyfm(x) as follows 



(TV) l 

V iN) = f- x \n(e^ (N) {x)) (3? 



we obtain for it a linear Schrodinger equation, i.e. 

d 2 N — 7 x 2 



& + — + t)W> = - (39) 

For even values of N we have the following simple solutions for this equa- 
tion: 







.2 

= , 


(40) 




;*) 


= e * , 


(41) 




;*) 


.2 

= xe 4 , 


(42) 


V (JV=2) 




= (l-X 2 )e"T, 


(43) 


^ (JV=0) 




/ 9 \ x ^ 

= x (3 — a; ) e 4 . 


(44) 



The above solutions for iV = 4, 2, are physical and therefore it is possible to 
calculate the following partial wave /i, for these cases 

I («> 

1/2 + -^, (48) 



,(JV= 


=4) 


a; 




,(/V= 


=2) 






,(iV= 


=0) 





+ vfe a; — Vb 



13 13 13 

-- — + '-^ + '—=. 47 

w w + V36 w- V36 

These poles in the w-plane lead to the following double-logarithmic asymptotic 
behavior of the corresponding scattering amplitudes: 

I 2 l<3| 2 

A UN) (s,t) = ^4° m (^) 9 ^ rW(s,t), (48) 

r 8+ioo j / \ w f (iV) 

r<">M) = / ^ ^) J —> (49) 



-lOO 
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where 



~ {N=i) (s,t) = 1, (50) 

^ 2) (M) = U(t^) h (51) 



2 \ Vkl 2 / Vk 12 

r(>{s ' t) = l[ 1+ (w) + (w) (52) 

The second solution of the Schrodinger equation for even N can be con- 
structed with the use of the Wronskian for two independent solutions 

*$)^ * % - *(3)^ *w = constant • ( 53 ) 

By integrating it we obtain the physical solution for M = 8 supergravity: 

f°° _al f°° 

ty(N=8)(x) = e 4 / &e 2 = e * / dye ? e xy , (54) 
./x Jo 

which is proportional to the probability integral erfc(\/2x). We can then obtain 
the corresponding scattering amplitude in the form 

ra+ioo j / 

r®(8,t) = - ^ A dye-Ve—. (55) 

Ja-ioo 27T2 \|g|V dx J 



Generally, the physical solution of the Schrodinger Eq. (39) can be ex- 
pressed in terms of the parabolic cylinder function 

e~T r°° dy _y* 

L y 



y { N)(x) = D^(x), D u (x) = —- -^c-Te"*". (56) 



Therefore 

fi N) 2 Id,/ r dy _,2 \ u 



a 6-NVbdx ln \l 7¥ e '' e ")' X= ^' <5?) 
The integral over y is convergent for v < 0. For example, for N = 7 we have 

/- 7) 2 d f r 00 dy _v2_ \ u 



_ } In / -^ e -^ e -** , ar = — . (5* 
w y/bdx \J y/y J y/b 
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In the general case with arbitrary v we can choose the integration contour 
L in the complex plane y which goes from y — +00, surrounds the singular 
point y — and returns again to y = +00: 



ti N) 2 1 d I 1 f dy _y* \ u 




co 6-N^bdx \2mJ L (- y )^ /' " y/b ' 

In particular, for N — > 6 we obtain 



(59) 



^ Vb Jo w 



00 



A6) 1 i-oo 2 1 

' dye-?e- xv = -\ dze- z e~t^, (60) 



which is in agreement with Eq. (34). After differentiating over a? in Eq. (59 ) and 
taking N = 4, 2, we can also reproduce our previous results in Eqs. ( |45[46|47fl . 

For odd values of N and N = 8 the function ^^(x) has an infinite num- 
ber of zeros situated asymptotically close to the lines argz = ±|7r. The 
trajectories of these Regge poles satisfy the following equation at large n 

X W 2 ~ _ 2 (7-N) lnx (A ° + 2 7re ± i^(2n + 1) , n = 0,1,2,.... (61) 

For M = 1 SUGRA there are three Regge poles with the real trajectories 

x[ N=1) » -2.460 , xf =1) « -0.452 , xf =1) = 1, 402 , (62) 

leading to the following growth of the scattering amplitude: 

9 / _ \ 1,402 ^\q\ 

limrM( S) ^ • (63) 

s^oo 5 \\q\ z J 

In Af = 3 SUGRA there are two real Regge trajectories, 

x[ N=3) « -1.747 , x^ =3) « 0.5508 , (64) 

also generating a growing contribution to the amplitude: 

2 / n 0,5508 y/a\q\ 

lim r (7V=3) (s,t) . (65) 

s-s>oo 3 \\q\ 1 ) 

For the case N = 5 there only exists one Regge pole with the real Regge 
trajectory 

X [ N =V w _o.762 , (66) 
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Figure 1: Scattering amplitude for N = 0,1, 2, 3, 4, 5, 6, 7, 8. 

which leads to an amplitude falling with energy (together with some oscilating 
contributions from the poles in the complex plane) 

i x -0.762 y/*\n\ 

lim r^ N=5 \s,t) «2 ^ . (67) 



Q 



In the cases N = 7,8 all the poles of f u are situated in the left hand side of 
the complex w-plane and therefore r^ N \s,t) here tends to zero when s — > oo. 
In particular, 

/ s \ -1-916 VS|g| / \ 
lim r {N=8) (s, t) = 2 — cos 2.8164 y/^\q\ In — . (68) 

s^oo \\q\ 2 J V \i\ J 

Besides these asymptotic estimates, we have performed an exact numerical 
analysis of the function r^(s,t) (see Eq. Q). For N = 0, 1, 2, 3, 4, 5, 6, 7, 8, 
in agreement with the previous analysis in this section, we show the energy 
behaviour of the scattering amplitudes for different supergravities in Fig. 0. 
The monotonically growing with energy solutions for N < 4 are shown in more 
detail in Fig. |2j The critical solution at iV = 4 is flat with energy The two 
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6 8 10 12 14 16 18 20 

s/(-t) 

Figure 2: Scattering amplitude for N = 0,1, 2, 3. 

monotonically decreasing with energy solutions for N = 5,6 are given in Fig. [3] 
and the two oscillatory and decreasing with energy ones for N = 7,8 can be 
seen in Fig. |4j 

It is interesting to point out that the falling asymptotic behavior of the 
amplitudes at s — > oo in the A^ = 5, 6, 7, 8 supergravities is likely to be re- 
lated to good ultraviolet properties of these theories, including the possible 
renormalizability of M = 8 SUGRA. 

3 Conclusions 

We have calculated the leading double-logarithmic in energy contributions to 
four-graviton scattering to all orders in the gravitational coupling. These terms 
are subleading with respect to eikonal contributions but important to under- 
stand the high energy behaviour of the scattering amplitudes. Our results are 
valid for any supergravity as well as for Einstein-Hilbert gravity. We have 
used infrared evolution equations which take into account both ladder and 
non-ladder topologies. The truncation of our resummation to two loops is in 
exact agreement with recent calculations in the literature for A/" = 4, 5, 6, 8 su- 
pergravities. Our results show a growth with energy for the amplitudes when 
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Figure 3: Scattering amplitude for N = 5, 6. 




Figure 4: Scattering amplitude for iV = 7, 8. 
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Af < 4, a critical invariance with the energy for Af = 4, and an asymptotic 
approach to zero when Af > 4. 
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